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CHAPTER  1 
INTRODUCTION 


The  ultimate  objective  of  this  work  is  to  obtain  accurate  predictions  of  the 
effect  of  structural  damage  on  the  collapse  strength  of  submarine  pressure 
hulls.  Therefore,  it  is  necessary  that  a  procedure  be  established  for  deter¬ 
mining  the  critical  static  buckling  load  of  typical  hull  structures  with  and 
without  pre-damage. 

Reference  1  used  the  finite  element  method  to  address  itself,  albeit  in  a 
narrow  way,  to  the  question  of  critical  static  collapse  pressure.  However, 
several  advanced  issues  were  not  treated  there.  These  relate  to  the  effect  of 
circumferential  boundary  conditions,  the  use  of  half  and  quarter  models,  the 
effect  of  mesh  size  on  eigenvalues,  plate  element  types,  the  effect  of  a  hole, 
and  finally,  the  effect  on  the  result  of  deviations  from  a  purely  cylindrical 
geometry.  It  is  the  intention  of  this  report  to  address  these  topics,  particu¬ 
larly  in  the  context  of  a  finite  element  buckling  analysis  which  presents 
considerable  practical  difficulties. 


Moussouros,  M. ,  Comparisons  of  Static  Collapse  Predictions  of  a  Ring 
Stiffened  Cylindrical  Shell  Subject  to  Hydrostatic  Pressure,  NSWC  TR  81-325, 
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CHAPTER  2 

BACKGROUND  TO  CRITICAL  LOAD  ANALYSIS 


In  the  following  paragraphs  we  briefly  discuss  the  phenomenon  of  buckling  of 
a  structure. 

When  a  structure  is  subjected  to  compressive  external  loads  there  always 
exists  a  critical  load  above  which  deformations  and  strains  increase  so 
drastically  that  the  load  carrying  capacity  of  the  structure  is  impaired.  This 
notion  of  the  critical  load  is  very  useful  to  damage  prediction. 

To  describe  this  instability  phenomenon,  one  may  form  the  potential  energy 
functional  of  the  discrete  equivalent  system  in  terms  of  the  generalized 
coordinates  (displacements)  and  a  loading  facto r  A  in  the  form 

V(q.,  A)  =  U(qi)  -  A  f(qi) 

where  U(q^)  is  the  strain  energy  of  the  structure  and  f(q^)  is  a  function  of 
the  deformations  (generalized  coordinates)  and  A  f(q^)  is  a  term  representing 
the  work  done  on  the  structure.  The  changes  of  such  a  functional  are  excellent 
indicators  of  the  behavior  of  the  system.  Moreover,  the  separate  contributions 
to  strain  energy  may  offer  further  insight  as  to  which  terms  dominate  as  the 
point  of  instability  is  approached. 

Geometrically,  the  maxima  and  minima  of  the  curve  relating  q^  and  A  are 
limit  points.  There  the  stability  of  an  equilibrium  path  changes  and  the  total 
potential  energy  has  an  inflexion.  In  general,  there  may  be  several  equilibrium 
paths,  with  intersections  referred  to  as  "bifurcation  points."  For  example,  in 
the  case  of  the  Euler  column  subject  to  an  axial  compressive  load,  the  primary 
path  and  the  secondary  path  cross  at  a  bifurcation  point  known  as  the  Euler 
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buckling  load.  Bifurcation  points  may  be  further  classified  as  either 
asymmetric,  stable -symmetric ,  or  unstable-symmetric,  depending  on  whether  the 
potential  energy  has  an  inflexion,  a  minimum  (second  variation  positive),  or  a 
maximum,  (second  variation  negative).  With  this  background  we  proceed  to  the 
modeling  process. 
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CHAPTER  3 

STRUCTURAL  MODELING 


We  starC  by  describing  the  steps  involved  in  an  analysis  by  the  finite 
element  method. 

(1)  Develop  a  model  of  the  real  structure,  say  an  internally  stiffened 
-cylinder,  by  subdividing  it  into  elements  (Figure  1)  of  various  shapes.  In 
attempting  to  model  the  expected  behavior  of  the  structure,  it  should  be  noted 
that  different  element  types  and  their  arrangement  may  alter  the  results 
drastically.  A  way  to  include  imperfections  must  be  devised.  This  may  be  a 
more  important  factor  in  buckling  computations  than  in  many  other  types  of 
struct'  al  analysis.  In  deciding  how  much  of  the  cylindrical  structure  to 
include  in  the  model,  we  note  that  we  could  consider  an  arbitrary  number  of  bays 
away  from  the  ends,  provided  they  are  equidistant.  We  may  impose  continuity 
conditions  at  arbitrary  cuts  for  symmetrical  structures  with  symmetrical 
loading.  In  the  present  instance,  it  appears  reasonable  to  use  a  minimum  of 
three  bays. 

(2)  Obtain  required  "inputs"  including  material  properties,  plate 
thicknesses,  properties  of  members,  etc. 

(3)  Simulate,  as  exactly  as  possible,  the  boundary  conditions,  whether  they 
are  of  the  stress  or  displacement  type.  Furthermore,  using  a  half  or  a  quarter 
model  of  the  structure,  the  computer  requirements  for  the  analysis  can  be 
reduced  considerably.  This  action,  however,  requires  additional  boundary 
conditions  at  the  interfaces,  and  these  may,  themselves,  create  difficulties 
which  must  be  considered  in  buckling  problems. 

(4)  Impose  the  external  loading  known  as  the  baseload.  This  can  then  be 
magnified  by  a  loading  factor.  Applying  the  external,  longitudinal, 
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compressive  forces  can  be  achieved  in  a  variety  of  ways.  One  may  include  heavy 
end  plates  which  model  the  end  closures  and,  simultaneously,  attempt  to  avoid 
their  buckling  and  then  enforce  the  pressure  load.  This  method  would  result  in 
heavy  local  bending  at  the  junctions,  possibly  affecting  the  buckling  mode. 
Another  method,  the  one  employed  in  this  report,  is  to  introduce  the  end 
pressure  as  a  uniform  axial  line  load  (Reference  2). 

(5)  Reconsider  items  (1)  -  (4)  when  the  shell  contains  a  hole. 


Moussouros,  M.,  Comparisons  of  Static  Collapse  Predictions  of  a  Ring  Stiffened 
Cylindrical  Shell  Subject  to  Hydrostatic  Pressure. 
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CHAPTER  4 
APPLICATION 


In  Reference  3,  a  portion  of  a  WW  II  submarine  was  analyzed  using  the 
4 

STAGS  general  shell  finite  element  program.  The  bifurcation  (linear  stress 
state)  option,  employing  the  full  length  of  three  frame  spaces,  but  only  a 
semicircular  portion  of  the  modeled  cylinder  was  used.  At  the  resulting  axial 
cuts  e  ■  0  and  6  “  ir  (Figure  2)  symmetry  conditions  were  employed.  It  was 
later  realized  that  such  an  analysis  might  not  yield  the  lowest  critical 
pressure. 

It  will  shortly  be  evident  that  mesh  size,  choice  of  plate  element,  and 
deviations  from  a  perfect  cylindrical  geometry  may  be  significant.  The  choice 
of  boundary  conditions,  such  as  the  ones  imposed  to  reduce  the  dimensionality  of 
the  analysis,  may  also  affect  the  results.  These  points  were  not  addressed  in 
Reference  5.  Therefore,  in  order  to  be  more  comprehensive,  irrespective  of  the 
presence  or  absence  of  predamage,  this  study  investigates: 

(1)  The  effect  of  boundary  conditions  used  to  reduce  dimensionality  of  the 
problem; 

(2)  The  effect  of  mesh  refinement  on  computed  eigenvalues; 


ntfoussouros,  M.,  Comparisons  of  Static  Collapse  Predictions  of  a  Ring  Stiffened 
Cylindrical  Shell. 

a 

Almroth,  B.  0.,  and  Brogan,  F.  A.,  and  Stanley,  G.  M. ,  Structural  Analysis  of 
General  Shells,  Vol.  II,  User  Instruction  for  STAGSC,  LMSC-D633873,  Apr  1979. 

^Moussouros,  M.,  Comparisons  of  Static  Collapse  Predictions  of  a  Ring 
Stiffened  Cylindrical  Shell. 
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(3)  Solutions  by  two  alternative  elements; 

(4)  Degradation  of  strength  due  to  a  hole;  and 

(5)  The  effect  of  initial  imperfections  on  the  computed  critical  pressure 
Comparisons  with  other  sources  and  programs  such  as  NASTRAN,  are  also  included 

1.  EFFECT  OF  BOUNDARY  CONDITIONS  IN  THE  MODELING  PROCESS 


Ue  establish  first  two  frames  of  reference,  a  global  Cartesian  (X,Y,Z)  and 
local  cylindrical  (r,8,x).  For  convenience,  the  boundary  conditions  will  be 
referred  to  the  local  system.  Furthermore,  we  define  the  following  terms 
pertaining  to  the  boundary  conditions: 


(1)  Type  I  (Symmetry): 


U 

U 


ee 


rr 


0 

0 

0 


(2)  Type  II  (Continuity): 


ee 


rr 


0 

0 


(3)  Type  III  (Symmetry): 


U. 


xx 


rr 


0 

0 

0 


(4)  Type  IV  (Antisymmetry): 


U 


0  (No  longitudinal  translation) 
0  (No  radial  translation) 

0  (No  rotation  about  tangential 
axis) 


(5) 


Type  V  (Relaxed  Antisymmetry): 


0 

0 


(6)  Type  VI  (Horizontal  Rigid  Body  Mode):  At  both  x"0  and  x*3s  and  e**!* 
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where  s  represents  to  frame  spacing  in  the  model 

(7)  Type  VII  (Horizontal  Rigid  Body  Mode):  At  both  x=0  and  x=3s  and 
0=3*/2 

U9  -  ° 

(8)  Type  VIII  (Vertical  Rigid  Body  Mode):  At  both  x=0  and  x=3s  and  0=0  and 
0=r 

ue  =  0 

In  order  to  investigate  the  effect  of  boundary  conditions  on  the  computation 
of  the  critical  pressure,  the  following  analyses  have  been  carried  out  on  a 
portion  of  a  circular  cylindrical  shell  encompassing  three  frame  spaces  (3s)  and 
four  equal  T-stiffeners  (rings):  These  models,  referred  to  as  models  (1)  -  (5), 
are  defined  as  follows: 

(1)  Full  cylinder  axially  and  peripherally  with  Type  I  boundary  conditions 
on  the  first  cross  section  and  Type  II  on  the  second.  To  prevent  vertical  and 
horizontal  rigid  body  motions,  we  apply  boundary  conditions  of  Type  VI,  VII,  and 
VIII.  (Boundary  conditions  (a)  of  Table  1.) 

(2)  Full  cylinder  axially,  but  only  one  half  circumferentially  with  Type  I 
boundary  conditions  on  the  first  cross  section  and  Type  II  on  the  second.  At 
the  generators,  or  cuts,  0  =  0  and  0  =  ir ,  apply  Type  III  boundary  conditions. 

To  remove  horizontal  rigid  body  motion  we  use  Type  VI  boundary  conditions  too. 
(Boundary  conditions  (b)  of  Table  1.) 

(3)  Full  cylinder  axially,  but  only  one  half  peripherally  with  Type  I  and 
Type  II  boundary  conditions  at  the  two  cross  sections  and  Type  III  and  Type  IV 
at  0  »  0  and  0  ■  ir.  Also,  apply  Type  VI  conditions.  (Boundary  conditions 

(c)  of  Table  1.) 

(4)  Full  cylinder  axially,  but  only  one  half  peripherally.  Instead  of  Type 
IV  conditions  at  0  ■  it,  as  in  Case  (c),  apply  Type  V  conditions.  (Boundary 
conditions  (d)  of  Table  1.) 
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(5)  Case  (2)  except  that  Type  III  conditions  at  the  cuts  6=0  and  6  = 
ir  are  replaced  by  Type  V.  (Boundary  conditions  (e)  of  Table  1.)  These  are 
modified  antisymmetry  boundary  conditions. 

а.  Description  of  Models 

Model  (1)  was  discretized  using  the  STAGS  410  quadrilateral  plate  element 
with  6  D.0.F.*/node.  The  inplane  translations  normal  to  the  edge  were  assumed 
linear,  while  the  translations  parallel  to  the  edge  were  cubic.  The 
out-of-plane  deformation  was  taken  as  cubic.  The  first  and  third  bays  were 
discretized  using  3  x  25  points,  where  3  indicates  the  number  of  nodes  axially, 
and  25,  the  number  of  nodes  peripherally  with  0  =  15°,  and  an  aspect  ratio  B, 
where 

B  ,  2R  -3'"{-9-/2>  =  9.0045/5  *  1.8 

and  R  is  the  shell  radius,  (34.50  in).  Ax  the  element  spacing  axially  (half  a 
frame  spacing)  and  e  “  15°.  The  middle  bay  has  a  higher  aspect  ratio  with 
Ax  *  2.5".  Hence  B  *  3.6. 

Models  (2)  -  (5)  were  analyzed  using  four  axial  subdivisions  on  the  first 
and  third  bays  and  eight  in  the  middle  one.  Circumferentially  there  were  13 
points  with  the  same  aspect  ratios  and  element  type.  In  addition  to  the 
previous  five  cases,  six  additional  combinations  were  also  used.  The  last  case 
(Case  11)  is  a  NASTRAN  computation  with  boundary  conditions  (b),  similar  to  Case 

б.  Computed  critical  pressures  are  given  on  Table  1. 

It  is  apparent  that  introducing  antisymmetry  at  6  *  ir,  while  maintaining 

symmetry  at  0  *  0  (Case  3)  (which  incorporates  the  boundary  condition  U  *  0 

2X 

at  0  ■  s),  lowers  the  critical  pressure  considerably  from  292.57  Lb/in  to 
o 

198.27  Lb/in  .  This  may  result  from  the  introduction  of  axial  reaction  forces 
at  the  antisymmetric  side.  This  is  further  supported  by  the  fact  that  for 


*D.0.F.  »  Degrees  of  Fteedom 
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Case  4,  which  omits  the  boundary  condition  U  *  0  at  e  *  it,  for  all  x  (0 

X  .2 

<  x  <  3s)  the  critical  pressure  is  increased  to  259.8  Lb/in  from  198.27 
2 

Lb/in  .  Similar  results  follow  for  Case  5  with  relaxed  antisymmetry  at  both 
ends  e  a  0  and  e  *  r. 

The  full  versus  half  model  comparisons  (Cases  1  and  2,  respectively) 

2 

indicate  a  small  reduction  in  the  critical  pressure  from  292.57  Lb/in  (Case  2 

2 

and  element  type  411)  (see  also  Reference  1)  or  295.50  Lb/in  (Case  6  and 
element  410)  to  263.38  Lb/in^. 

Table  1  also  includes  results  for  other  combinations.  Case  9,  for  example, 

differs  from  Case  7  in  the  aspect  ratio.  Mesh  refinement  to  a  more  balanced 

2 

aspect  ratio  results  in  bypassing  the  lowest  eigenvalue,  286.05  Lb/in  ,  and 

2 

converges  instead  to  440.71  Lb/in  .  This  is  discussed  later.  The 

corresponding  mode  shapes  for  Cases  1-5  are  given  in  Figures  2  through  11.  As 

has  already  been  mentioned,  the  model  of  Case  2  or  6  (which  only  differ  in  the 

plate  element  type)  and  is  reported  as  Case  11  is  run  using  NASTRAN.  The 

2 

critical  pressure  is  obtained  as  390.46  Lb/in  .  It  is  close  to  the  circular 

2 

ring  solution  reported  in  Table  2  as  391.63  Lb/inS 

b.  Effect  of  Initial  Pressure 

2 

For  all  the  above  cases,  the  applied  baseload  pressure  was  1  Lb/in  .  Case 

2 

1  was  then  recomputed  using  a  baseload  pressure  of  100  Lb/in  .  The  resulting 
critical  pressure  and  modes  were  identical,  suggesting  independence  of  the 
result  from  this  aspect  of  the  procedure. 

c.  Effect  of  Live  Load 

For  all  cases  reported  here,  unless  otherwise  noted,  the  pressure  load  was  a 
"dead  load,"  i.e.,  did  not  change  direction  of  action.  Case  1  (full  model)  was 
run  with  STAGS'  "live  load"  option  for  pressure,  yielding  4p  =  280.14 

2  2  Cr 

Lb/in  compared  to  the  dead  load  result  of  263.38  Lb/in  .  The  corresponding 
critical  mode  was  identical  to  that  of  Figure  2.  Appendix  A  contains  Tables  A-l 
and  A-2  displaying  the  mode  1  displacement  distributions  for  the  dead  and  live 
pressures  of  Case  1. 


13 
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FIGURE  4. 


MODE  1  FOR  CASE  2,  HALF  MODEL  PERIPHERALLY,  SYMMETRICAL  BOUNDARY 
CONDITIONS,  Apcr  -  292,57  LB/IN2 
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FIGURE  5.  PLOT  OF  ORIGINAL  FORM  AND  MODE  1  AT  A 


NSWC  TR  82-172 


FIGURE  7.  PLOT  OF  ORIGINAL  FORM  AND  MODEL  1  AT  A  CROSS  SECTION  AT  x  -  1 .5s,  CASE  3, 
VERTICAL  SCALE  IS  TWICE  HORIZONTAL,  MAGNIFICATION  -  4 


18 
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d.  Mode  Shapes 

Examination  of  the  mode  shape  corresponding  to  the  calculated  critical 
pressures  suggests  that  the  ring  reinforced  shell  will  buckle  in  an  asymmetric, 
local,  mode  otherwise  referred  to  as  lobar  or  diamond-shape  buckling.  The  known 
forms  are:  (1)  Axisymmetric  local  buckling  (combination  of  yielding  and 
buckling),  (2)  Asymmetric  local  buckling  in  between  the  ring  stiffeners,  and  (3) 
General  Asymmetric  buckling.  These  forms  are  displayed  in  Figures  12  through 
14.*  They  correspond  to  bifurcation  rather  than  limit  points,  as  already 
mentioned. 

e.  Comparison  with  Other  Computations 

Table  2  is  a  summary  of  results  obtained  with  a  computer  program  based  on 
various  formulae  for  elastic  buckling.  We  note  that  "Fluegge's  formula"  *  is 
for  general  instability,  while  "Fluegge's  formula  (skin  only)"  is  for  interbay 
instability.  Most  of  these  results  were  reported  in  Reference  8,  where  they 
were  calculated  manually.  These  can  also  be  compared  with  the  NASTRAN  result 
reported  in  Table  1  and  discussed  previously. 


♦Reproduced  from  Renzi,  J.  R.,  Optimization  of  Orthotropic,  Non-Linear,  Ring- 
Stiffened  Cylindrical  Shells  Under  External  Hydrostatic  Pressure  as  Applied 
to  MMC  Materials,  NSWC  TR  79-305,  30  Sep  1979. 

^Fluegge,  W.,  "Die  Stabilitaet  der  Kreiszylinderschale,"  Ingenieur  Archiv, 

Vol.  3,  pp.  403-506,  1932. 

^Timoshenko,  S.  P.,  and  Gere,  J.  M.,  Theory  of  Elastic  Stability  (McGraw-Hill 
Book  Co.,  Second  Edition,  1961),  pp.  457-519. 

^Moussouros,  M.,  Comparisons  of  Static  Collapse  Predictions  of  a  Ring  Stiffened 
Cylindrical  Shell. 
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FIGURE  12.  AXISYMMETRIC  LOCAL  BUCKLING 


& 
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FIGURE  14.  GENERAL  ASYMMETRIC  BUCKLING 
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TABLE  2. 

SUMMARY  OF  CRITICAL  PRESSURES  AS  PER  DIFFERENT 
METHODS  IN  THE  OPEN  LETERATURE 


CIRCULAR  RING  FORMULA.  CRITICAL  PRESSURE  =  391.639  LB/SQ.IN. 
R.  VON  MISES  FORMULA.  CRITICAL  PRESSURE  =  494.266  LB/SQ.IN. 
SOUTHWELLS  FORMULA.  CRITICAL  PRESSURE  =  380.569  LB/SQ.IN. 
U.S.  EXP.  MOD.  BASIN  FORMULA.  CRITICAL  PRESSURE  =  495.622  LB/SQ.IN. 
FLUEGGES  FORMULA.  CRITICAL  PRESSURE  =  2161.224  LB/SQ.IN. 
FLUEGGES  FORMULA  (SKIN  ONLY.)  CRITICAL  PRESSURE  =  494.230  LB/SQ.IN. 


NOTE:  Inertia  of  plate  stiffener  combination  is  calculated  assuming  full 
effectiveness,  I  =  1.7869  in4. 
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f.  Stress  State  After'  Buckling 

Examination  of  the  resulting  stresses  corresponding  to  a  critical  pressure 
2 

AP  =  263.68  Lb/in  gives  at  midbay  2  a  longitudinal  stress  a  = 
rcr  °  xx 

-24.15  ksi  (i.e.  47%  of  yield),  a  circumferential  stress  o9Q  =  -37.92  ksi 
(75%  of  yield)  (a„v/oaa  =  0.637),  and  a  maximum  stiffener  stress, 

-20.62  ksi  (40%  of  yield).  The  buckled  configuration  exhibited  only  elastic 
deformation. 

g.  Comments  on  Additional  Computations 

Additional  computations  have  been  performed  using  STAGSC  to  evaluate  a 

q 

method  proposed  by  Brendel  and  Ramm.  According  to  these  authors,  one  first 
obtains  the  approximate  number  of  circumferential  waves  (m)  through  a  formula 
such  as 


or  by  using  Figure  15.  In  the  present  case,  the  formula  gives  m  =  17.54,  while 
Figure  15  indicates  either  m  =  14  or  m  *  15.  For  the  model  results  reported 
below,  the  effects  of  ring  stiffeners  were  incorporated  using  discrete  beam 
elements . 

These  values  are  only  approximate  and  it  is,  therefore,  necessary  to  compute 
the  minimum  eigenvalue  for  a  variety  of  angles,  corresponding  to  m  peripheral 
waves,  and  for  suitable  aspect  ratios.  Table  3,  formed  by  computing  6 
segments  from  18°  to  180°  (m  =  20  to  m  =  2),  indicates  that  m  =  3  corresponus.  to 
the  critical  pressure,  the  magnitude  of  which  agrees  fairly  well  with  the  full 
and  even  the  half  models  (Cases  1  and  2).  There  is,  however,  disagreement 
between  the  number  of  waves  at  buckling.  Figure  3  and  Table  A-l  tend  to  support 


9 

Brendel,  B.,  and  Ramm,  E.,  "Linear  and  Nonlinear  Stability  Analysis  of 
Cylindrical  Shells,"  Journal  of  Computers  &  Structures,  Vol.  12,  pp.  549-558, 
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TABLE  3.  BUCKLING  ANALYSIS  OF  PANELS  BETWEEN  RING  STIFFENERS 


6(0) 

[One  Wave] 

No.  of 

Segments 

Peripherally 

No.  of 
Segments 
Axially 

Aspect 

Ratio 

As/Ax 

/min  ? 
(Lb/in^) 

180 

12 

4 

0.900 

308.4 

120 

12 

2 

1.203 

259.0 

90 

6 

2 

1.800 

484.1 

60 

6 

4 

2.406 

561.7 

24 

6 

4 

1.440 

782.6 

22.5 

6 

4 

1.351 

773.9 

18 

6 

4 

1.083 

813.1 
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6  half  waves  over  Che  entire  ring  (2  half  waves  over  120°)  versus  1  half  wave  by 
the  previous  process.  See  Table  A-3.  Figure  16,  which  shows  the  mode  1  shape 
for  a  panel  with  an  arc  8  3  120°  is  very  similar  to  that  of  Figure  2.  Figure 
17  shows  mode  1  for  the  case  of  a  panel  e  =  90°,  assumed  to  contain  1  full 
wave.  The  computation  indicates  5  half  waves  over  the  90°  arc.  On  the  face 
of  it,  the  previous  method  may  be  unreliable  as  far  as  mode  shape  is  concerned. 
It  must  be  emphasized,  however,  that  aspect  ratio  has  a  marked  effect  on  the 
critical  load  and  the  corresponding  mode  as  is  evidenced  in  Table  4  for  the 
panel  analysis. 

For  most  cases,  except  8  =  180°  and  120°,  as  the  aspect  ratio  goes  up 
(number  of  circumferential  elements  versus  axial  ones),  Apc^  is  reduced.  Note 
that,  if  the  aspect  ratio  is  reduced  considerably,  the  critical  pressure  in¬ 
creases  substantially  in  the  direction  of  that  of  a  simply  supported  strut  of 
length  s*  between  ring  stiffeners  subject  to  axial  loading. 

2.  EFFECT  OF  MESH  REFINEMENT  ON  EIGENVALUES 

The  following  discussion,  though  not  comprehensive,  appears  to  be  generally 

»* 

applicable.  (See  Table  5.)  That  is,  a  refinement  of  the  mesh  in  the  circum¬ 
ferential  direction  (increase  of  aspect  ratio  from  3.60  to  7.20)  lowers  the 
computed  eigenvalue  from  295.00  Lb/in  to  286.86  Lb/in  (Case  1  to  Case  2). 
Refinement  by  increasing  the  number  of  elements  by  a  factor  of  two  (Case  1  to 
Case  3),  while  maintaining  essentially  the  same  aspect  ratio,  leads  to  confusing 
results.  It  appears  that  mesh  refinement  tends  to  increase  the  dimensionality 
of  the  problem  to  the  point  at  which  accurate  extraction  of  the  lowest 
eigenvalue  becomes  difficult. 


2  2  2  6  2  2 
*o  »  *  E (il)  -  *  x  30  x  10  x  (Oj.2A5)  3  14,810  Lb/in 

CR  T7  s  IT  10 
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TABLE  4. 

BUCKLING  ANALYSIS  OF  PANELS  BETWEEN  RING 
STIFFENERS  FOR  A  RANGE  OF  ASPECT  RATIOS 


No.  of 

No.  of 

As 

No.  of  Waves 

0(°) 

Segments 

Segments 

Ra 

Peripherally 

[one  wave] 

Peripherally 

Axially 

AS 

2  3( 

2  3( 

)  6  2  1 

)  12  2  1 

5  2  i 

5  2i 

\  6  4 

\  12  2  1 

6  22.5  6  4 

6  .  22.5  12  2  1 

18 

6 

4  1 

18 

12 

2  C 

(lb/in2) 


561.7 
1 ,116.5 
1,089.5 


2,862.9 

2,965.0 


773.9 

2,650.3 
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3.  SOLUTIONS  BY  ALTERNATIVE  ELEMENTS 

The  present  version  of  STAGSC  offers  a  choice  between  plate  elements  series 

410,  411,  412  and  420,  421,  422.  Since  we  have  already  used  410,  we  employed 

411  as  a  comparison.  The  quadrilateral  plate  element  411  has  32  D.O.F.  with  3 

translations,  2  inplane  rotations,  2  independent  normal  rotations  at  each 

corner,  and  1  tangential  displacement  at  each  mid-side  node.*®  From  Reference 

11  it  is  evident  that  there  are  difficulties  with  the  use  of  finite  elements  in 

instability  phenomena.  A  case  similar  to  Case  2  (Table  5),  for  the  aspect  ratio 

displayed  on  Table  1  as  Case  2,  and  employing  the  411  element  was  run.  The 

2 

results  using  these  different  elements  were  comparable  (295.00  Lb/in  versus 
2 

292.57  Lb/in  ).  On  the  face  of  it,  it  appears  that  for  the  type  of  problem 
considered,  element  410  is,  perhaps,  more  than  sufficient,  and  the  additional 
expense  of  the  extra  D.O.F.  of  element  411  is  not  warranted. 


*®Thomas,  K.,  and  Sobel,  L.  H.,  Evaluation  of  the  STAGSC  -  1  Shell  Analysis 
Computer  Program,  Westinghouse  Electric  Corporation,  Report  N.  Ward-10881,  Aug 
1981. 

**Almroth,  B.  0.,  and  Brogan,  F.  A.,  "Computational  Efficiency  of  Shell 
Elements,"  AMD-Vol.  48,  Nonlinear  Finite  Element  Analysis  of  Plates  and 
Shells ,  The  American  Society  of  Mechanical  Engineers,  Washington,  D.C.,  15-20 
Nov  1981.  *  | 
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4.  DEGRADATION  OF  STRENGTH  DUE  TO  A  HOLE 

References  12  and  13  include  some  results  of  structures  with  holes#  In  the 
present  analysis,  the  half  model  was  studied  with  a  cutout  on  one  side  only. 

The  hole  dimensions  were  ax*as.  The  actual  peripheral  dimension  in  the 
cutout  is  2As.  Computed  results,  employing  element  410,  indicate  that  there 
is  a  reduction  in  the  elastic  buckling  load  with  holing.  This  is  consistent 
with  results  of  Reference  14.  Table  6  summarizes  these  results.  Examination  of 
this  Table  suggests  the  danger  of  using  non-uniform  meshes  in  such  computations 
as  the  lowest  eigenvalue  may  be  missed  (Case  4).  Although  not  fully  conclusive, 
the  first  and  second  cases  indicate  a  reduction  in  strength  by  5%  because  of  the 
hole.  As  the  size  of  the  opening  is  increased  (from  Case  1  to  2  and  to  3)  the 
critical  pressure  is  reduced  correspondingly.  Figures  18-20  show  the  mode 
shapes  for  Case  2,  while  Figures  21  and  22  exhibit  the  critical  mode  shape  for 
Case  3.  Finally,  the  hole  has  a  mild  effect  on  the  stress  redistribution  with 
the  appearance  of  a  small  shear  stress  a  and  an  approximate  magnification 
of  longitudinal  stresses  by  1.5%  and  circumferential  ones  by  3.5%. 


5.  THE  EFFECTS  OF  IMPERFECTIONS  ON  CRITICAL  PRESSURE 


A  few  numerical  analyses  have  been  carried  out  with  the  imperfections  option 


of  STAGS  whereby  the  magnitude  of  the  amplitude  (Wq),  the  location  of  the 


maximum  amplitude  (X,d),  and  half  wave  lengths  of  the  imperfection  axially 
(X^)  and  tangentially  (eT)  were  taken  as 


^"^Brogan,  F.,  and  Almroth,  B.  0.,  "Buckling  of  Cylinders  with  Cutouts,"  AIAA 
Journal ,  Vol.  8,  No.  2,  pp.  236-240,  Feb  1970. 

^^Brogan,  F.  A.,  and  Almroth,  B.  0.,  "Practical  Methods  for  Elastic  Collapse 
Analysis  for  Shell  Structures,"  AIAA  Journal,  Vol.  9,  No.  12,  pp.  3421-2326, 
Dec  1971. 


^^Brogan,  F.,  and  Almroth,  B.  0.,  "Buckling  of  Cylinders  with  Cutouts.". 
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TABLE  6. 

EFFECTS  OF 

CUTOUTS  ON  CRITICAL 

PRESSURE 

Type  of  Model 

Element 

^Pe 

Half 

Opening  Size 

AX  *  AS 
(in  x  in) 

Same  as  Case  7,  (5  x  13), 
No  opening 


0x0 


Half  model  with  symmetry 
at  e  =  o,  0  *  ir.  Uniform 
mesh  in  middle  portion 
and  twice  as  fine  at  ends 
(5  x  13) 


2.5  x  9.032 


Same  as  Case  2  (3  x  13) 


5.0  x  9.032 


Half  model  with  symmetry 
at  0  =  o,  0  =  ir. 
Nonuniform  mesh  peri¬ 
pherally  and  axially 


0.50  x  0.25 


1 .25  IN,  A  $  -  9.0045  IN,  ( A  p)2  -  308.45  LB/IN2 


0 
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FIGURE  22.  STIFFENED  CYLINDRICAL  SHE  LL  WITH  OPENING.  (SAME  AS  FIGURE  21.) 
COMBINED  DISPLACEMENTS  (u,  v,  w)  MODE,  a  p  .  266.26  LB/IN2 
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=  O.lh  -  0.0245  in.,  XT 

O  L 


=  10  in.,  eL  =  =  10‘ 


and 


X  *  1.5s  ■  15  in.,  6  *  0  in  Cases  1,  2,  6  =  10°  in  Case  3. 


The  results  are  shown  in  Table  A-l.  Figures  23  through  27  display  the 
corresponding  critical  mode  for  the  case  analyzed. 


We  can  further  investigate  this  question  using  results  from  Budiansky  and 
Amazigo^  and  Bushnell^  relating  to  asymptotic  theory.  In  fact,  using 
Figure  28  for  the  corresponding  Batdorf  non-dimensional  parameter 


Z 


with  L  *  10  in,  R  *  34.5  in,  h  =  0.245  in,  v  =  0.300,  Z  =  11.285,  and  with 

hydrostatic  pressure  and  ring  stiffeners,  we  obtain  parameter  b  =  -  0.364. 

Since  W  *  O.lh  ■  0.0245  in.  the  Bushnell  formula,  where  AP  is  the 
o  s 

collapse  pressure  with  imperfections,  and  APC  of  the  ideal  structure,  is 


AP 
_ s 

AP 

c 


l-3(-  |> 


1/3 


P», 


and  assuming  p 


1,  we  obtain 


AP 

c 


0.886 


which  is  about  the  same  value  obtained  by  using  Figure  29.  This  results  in  a 

2  2 

critical  pressure  APg  *  261.8  Lb/in  ,  using  aP£  ■  295.5  Lb/in  without  the 


^Budiansky,  B.  E.,  and  Amazigo,  J.  C.,  "Initial  Post-buckling  Behavior  of 
Cylindrical  Shells  Under  External  Pressure,"  Journal  of  Mathematics  and 
Physics ,  Vol.  47,  pp.  223-235,  1968. 

^Bushnell,  D.,  "Buckling  of  Shells-Pitfalls  for  Designers,"  AIAA  Journal, 
Vol.  19,  No.  9,  pp.  1183-1226,  Sep  1981. 
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23.  FULL  STIFFENED  CYLINDRICAL  SHELL  WITH  IMPERFECTIONS,  (u,  v.  w)  COMBINED 
DEFORMATION  MODE.  AMPLITUDE  OF  IMPERFECTION  -  0. 1h.  MAXIMUM 
IMPERFECTION  AMPLITUDE  AT  x  -  IS*,  8  -  0,  Ap„  -  283.88  LB/IN2 
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FIGURE  24.  STIFFENED  CYLINDRICAL  SHELL.  SYMMETRY  CONDITIONS  AT  0-0.  0-  ir. 
(w)  RADIAL  DEFORMATION  MODE.  IMPERFECTION  AMPLITUDE  -  O.lh. 
MAXIMUM  IMPERFECTION  AT  x- 1.5i,  0-0,  A  p„  -  316.04  LB/IN2 
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FIGURE  26. 


STIFFENED  CYLINDRICAL  SHELL  WITH  IMPERFECTIONS.  SYMMETRY  CONDITIONS  AT 
8-0.  0-  JT  .  (w)  RADIAL  DEFORMATION  MODE.  IMPERFECTION  AMPLITUDE  -  0. 1  h 
MAXIMUM  IMPERFECTION  AMPLITUDE  AT  x  -  1.5*,  0- 10°,  Ap..  -  316.04  LB/IN2 
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imperfection.  This  compares  with  283.88  Lb/in  of  Table  7,  and  confirms  the 
expected  result  that  imperfections  lower  the  critical  pressure.  (Furthermore, 
it  suggests  that  263.38  Lb/in2,  as  obtained  by  the  full  model  (Table  1,  Case 
1),  is  an  artificial  characteristic  root.) 

6.  GENERAL  COMMENTS 

We  note  that  the  presence  of  the  hole  (Cases  2  and  3  of  Table  6)  appears 
to  remove  the  lowest  eigenvalue  obtained  in  the  model  without  a  hole  (see 
Table  1,  Case  1),  which  appears  to  be  a  spurious  result.  The  present  results 
(272.94  Lb/in  ,  266.26  Lb/in  )  compare  well  with  the  second  one  of  Case  1 
Table  1  (ap  ■  293.74  Lb/in2).  Similarly,  Cases  2  and  3  (Table  6)  compare 
with  Case  2,  Table  1,  in  which  the  imperfection  (Case  1,  Table  7)  again  seems 
to  have  caused  the  lowest  eigenvalue  to  be  missed. 
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TABLE  7. 

EFFECT  OF  IMPERFECTIONS  ON  CRITICAL  PRESSURE 


Type  of  Model 


Element 
T 


tj'ld 


Mesh  in 
Middle 
Bays 


Full  Cylinder  as  in 

Case  1,  Table  1,  with 
imperfections  (given 
as  Case  10  on  Table  1) 

410 

5  x  25 

283.88 

Half  Cylinder  as  in 

Case  2,‘  Table  1 ,  with 

410 

5  x  13 

316.04 

Imperfections  (0  *  0) 

V 

* 

As  Case  2  except  that 

angular  location  of 
maximum  imperfection 
at  0  *  10° 

410 

5  x  13 

316.04 
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CHAPTER  5 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  WORK 


In  this  report  only  the  bifurcation  option  (linear  stress  state)  of  STAGS 

was  studied  by  analyzing  a  portion  of  a  ring  stiffened  cylinder  subject  to 

2 

hydrostatic  loading.  The  critical  pressure  obtained  was  about  280  Lb/in  . 

The  following  conclusions  may  be  drawn  from  all  these  studies: 

(a)  Employ  the  analytic  formulae  to  obtain  an  approximate  value  for  the 
critical  load.  Start  at  a  higher  baseload.  Analyze  the  entire  cylinder 
peripherally,  but  half  axially.  If  such  a  model  is  too  big,  then  use  half  the 
model  peripherally  and  employ  symmetry  conditions  there.  Treat  the  resulting 
answers  as  approximate. 

(b)  Do  not  employ  a  very  fine  mesh,  as  it  may  tend  to  bypass  the  critical 
pressure. 

(c)  Element  410  is  as  good  as  Element  411  for  buckling. 

(d)  A  hole  appears  to  degrade  the  critical  load  slightly.  The  actual 
amount  depends  on  the  hole  size. 

(e)  If  imperfections  are  used,  it  is  preferable  to  use  the  full  cylinder 
model. 

In  the  use  of  circular  cylinders  subject  to  hydrostatic  loading,  both  the 
structure  and  the  load  are  symmetrical  and  the  prebuckling  behavior  is  non¬ 
linear.  Collapse  may  occur  in  an  axisymmetric  way  at  a  limit  point  or  in  an 
asymmetric  mode  if  bifurcation  occurs.  These  modes  can  be  analyzed  by  carrying 
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out  a  nonlinear  stepwise  analysis  under  increasing  load  (with  or  without 
imperfections).  Another  possibility  is  to  carry  out  a  dynamic  analysis  of  the 
same  structure  and  loading  over  a  time  much  larger  than  the  natural  period  of 
excitation  of  the  structure,  and  thereby  obtain  a  load  deflection  path.  For  the 
last  one,  estimates  of  the  natural  period  of  excitation  are  necessary. 
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TERMS 


Translation  radially  outwards 

Translation  in  longitudinal  direction  (X-axis) 

Translation  along  the  tangential  axis 

Rotation  about  the  radial  axis 

Rotation  about  the  longitudinal  axis  (X-axis) 

Rotation  about  the  tangential  axis 


TABLE  A-l.  MODE  1  TRANSLATIONS  (u,v,w)  FOR  DEAD  LOAD  PRESSURE 
CASE  1,  TABLE  1.  AT  x  =  1.5s,  Ap  =  263.38  Lb/in2 
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TABLE  A- 3.  MODE  1  TRANSLATIONS  (u.v.w)  FOR  DEAD  LOAD  PRESSURE  ON 
PANEL  OF  ARC  0°  -  120°  IN  THE  MIDDLE  OF  BAY  AS  PER 
TABLE  3,  FIGURE  17.  Ap  -  259.0  Lb/in2 


NSWC  TR  82-172 


OQOOOO^OO 

ooooooooo 

♦♦♦♦♦♦!♦♦ 

UitiJUlLklUJUJuJUJtiJ 

NaxoiMVlflS^-H 

co^ir>a»o'innm^ 

O'O'aJ'O'S-rvj— «(\jsf 


o  o  o  o 
o  o  o  o 
♦  ♦  ♦  ♦ 
UJ  U  UJ  LJ 
N  ^  ^  ^ 
(0>0  8  + 
in  — <  o»  m 
o  in  >♦  oo 

<0  35  i>  O' 


ooooooooooo 

OOOQOOOOOOO 
♦  ♦♦♦♦♦♦♦♦♦♦ 
UJliJUJUJUJUJtiJUJUJUiUJ 
OClHIflHSHlflHrlO 

xir^nsoanasma 

^fn^in^vo^inst-nr^ 


mmnnnrnmmmnn 

ooooooooooo 

•  •  l  l  i  l  i  i  i  i  i 

UJUJUJijJUJUiUiUJUJUJUi 

wommo'-d-.o— »;r>ror'. 

D  O  S  (\J  r-.  CT*  £\J  lO  O 

tO  40  ^  O'  lO  N  SO  ^  >*4  *— 4  | 
'C-4’  80'84’lfl-tfflO0> 
tfli04n(viHMM(u««i 


OOOOOOOOOOOOO 

OOOOOOOOOOOOO 

OOOOOOOOOOOOO 

•  •••••••••••a 

>-0000000000000 

noncNoo'OHivin^in 


OOOOOOOOOOOOO 

OOOOOOOOOOOOO 

OOOOOOOOOOOOO 

•  ••••••••••a* 

xintrirtintntntrttnintntniniD 


O  rm 4  _•  _| 

u 


(V  (\)  (\l  ru  (\l  fU  <\l  at  (M  fU  (M  (V  (\J 


NSWC  TR  82-172 


DISTRIBUTION 


Copies 

Chief  of  Naval  Material 
Attn:  CAPT  G.  Jarrett  (ASW-14)  1 

Dr.  A.  J.  Faulstich  (MAT  07)  1 

Navy  Department 
Washington,  D.C.  20360 


Copies 

Commander 

Naval  Weapons  Center 
Attn:  Code  533  (Technical 

Library)  1 

China  Lake,  California  93555 


Commander 

Naval  Sea  Systems  Command 
Attn:  SEA-03B 

SEA  63R  (F.  Romano) 
SEA-63R32  (Murphy) 
SEA-9G32 

SEA-32R  (C.  Pohler) 

SEA-3221  (H. .Ward) 

PMS-402 

PMS-406 

PMS-407 

Department  of  the  Navy 
Washington,  D.C.  20362 


Commander 

David  W.  Taylor  Naval  Ship 
Research  &  Development  Center 
Attn:  Code  17  (Dr.  W.  Murray) 
Code  175  (J.  Sykes) 

Code  175.2  (B.  Whang) 
Code  175.2  (T.  Giacofci) 
Code  175.2  (T.  Gilbert) 
Code  175.3  (W.  Conley) 
Code  175.3  (P.  Manny) 
Code  184.4  (M.  Hurwitz) 
Code  172  (M.  Krenzke) 

Code  1620.3  (R.  Jones) 
Betheada,  Maryland  20084 


1 

1 

1 

1 

1 

1 

1 

l 

1 


l 

1 

1 

1 

1 

1 

1 

1 

1 

1 


David  W.  Taylor  Naval  Ship  Research 
and  Development  Center 
Attn:  Code  177  (R.  Fuss) 

Code  177.1  (V.  Bloodgood) 

Code  177.1  (M.  Riley) 

Code  177. 1(R.  Higginbotham) 
Underwater  Explosion  Research  Division 
Portsmouth,  Virginia  23709 


Commander 

Naval  Ocean  Systems  Center 

Attn:  Technical  Library  1 

San  Diego,  California  92152 

Commanding  Officer 
Naval  Underwater  Systems  Center 
Attn:  (D.  J.  Lepore)  1 

Newport,  Rhode  Island  02840 

Office  of  Naval  Research 

Attn:  Code  474  (D.  N.  Basdekas)  1 

8t*0  North  Quincy  Street 
Arlington,  Virginia  22217 

Defense  Nuclear  Agency 

Attn:  SPSS  (LT  D.  Sobota)  1 

Washington,  D.C.  20305 

Defense  Technical  Information 
Center 

Cameron  Station 

Alexandria,  Virginia  22314  12 

Library  of  Congress 

Attn:  Gift  and  Exchange  Division  4 

Washington,  D.C.  20540 

Internal  Distribution: 


R10  1 
R102  2 
R14  20 
R14  (M.  Moussouros)  1 
E431  9 
E432  3 
E35  1 


(1) 

i..‘  m..»>  .  fc- ^-.v. 


49 


